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We have mapped the physics of a system of semi-fiexible inextensible polymers onto a complex 
Ginzburg-Landau field theory using techniques of functional integration. It is shown in the limit of 
low number density of monomers in a melt of semi-flexible, inextensible polymers, kept apart by an 
excluded volume interaction, the radius of gyration scales as A''"', where A'^ is the chain length, and 
v = 1, in contrast to the value of f « 0.6 for flexible polymer melts. Using Renormalization Group 
arguments, we show that the system exhibits an infra-red stable fixed point which can be identified 
as the transition to the entangled state. Experiments to test these calculations are suggested. 
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Polymers in a biological setting occur with a variety 
of persistence lengths, ranging from being much smaller 
than the overall length lif the polymer, to being compara- 
ble to the chain lengthtl. When the persistence length is 
much smaller than the chain length, the polymer can be 
described as being flexible. Generally however, the per- 
sistence length is not much smaller than the chain length, 
and the theoretical treatoients of such semi-flexible poly- 
mers becomes complexBfl, especially when one takes into 
account the fact that these polymers are inextensible, 
and can interact with each other. 

One purpose of this paper is to understand the physics 
of melts of inextensible, semi-flexible polymers in the low 
number density Jimit, by generalizing de Genne's theo- 
retical approaclD. In this paper, it will be shown that the 
Flory exponent v, which characterizes the radius of gyra- 
tion Rg ^ N'^ , where N is total number of monomers in 
an average chain, is different for the case of interacting 
semi- flexible polymers, and is given by z/ = 1. 

Renormalization Group arguments will also be used to 
show that as the number density of semi-flexible polymer 
segments increases in a melt, the system actually begins 
to behave as if the chains were becoming increasingly 
independent, analogous to Flory 's theorem for melts of 
flexible chains. Furthermore, it can be shown that the 
system scales towards a thermodynamic fixed point as 
the number density is increased further, indicating a sec- 
ond order phase transition to the entangled state. 

Treatments of semi-flexible polymers can become qwitf 
sophisticated, even when restricted to a single chainBLl. 
In fact the treatment of ensembles of interacting semi- 
flexible chains is an open problem. In this paper, the fo- 
cus will be on semi-flexible polymers interacting through 
hard-core repulsion between the monomers. This inter- 
action will be represented by an excluded volume term. 
It turns out that the elegant path integral formalism for 
a single semi-flexible chain put forth by KleinertQ is well- 
suited for generalization to a collection of interacting 
polymers. In this generalization, the physics of semi- 
flexible polymers gets mapped onto a time-dependent 
complex Ginzburg-Landau complex field theory. 

Time in the theory is merely a label for monomers along 
the chain. 



Let us begin with the worm-like model of Kratky and 
Porod for a single semi-flexible polymer: 
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where R{n) is the location of the n monomer in the 
chain which is N units long, t{n) is the tangent vector 
at R{n), (3 — 1/kT, k is Boltzmann's constant, T is the 
temperature, k = ip/lJ^, and £p is the persistence length 
over which the tangent vectors remain correlated, and h 
is the monomer length. Note that Eqn.(|^) has been writ- 
ten in a form that is reminiscent of that for the flexible 
chain case, with the replacement R t. 

Inextensibility can be incorporated formally as done 
by KleinertQ by introducing a correlation function via a 
path-integral representation. It expresses the fact that in 
(t, n) space, segments q units apart on the chain are dis- 
tributed in a Gaussian fashion (see Eqn.(^) below). This 
is esentially the definition of the persistence length along 
a semi-flexible chain: 
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where d — 3 refers to the dimension of the tangent vec- 
tor, and integration over the reduced (d — 1) dimensions 
incorporates a length preserving constraint. More specifi- 
cally, the restriction of the integration to d— 1 = 2 dimen- 
sions indicates the incorporation of a local inextensibility 
constraint, viz., i(n)^ — [dx / ds)"^ + {dy / ds)'^ + {dz / ds)'^ — 
constant. This restriction indicates that the integration 
is to be performed over the surface of a sphere in tangent 
space. This can lead to cumbersome calculations. An 
approximation will now be introduced, based on physical 
considerations, which greatly reduces the level of com- 
plexity of the subsequent calculations. 

One can imagine that if we have a collection of worm- 
like chains, such that it costs energy to bend any one 
of them, and their number density is quite low, it 
would be quite easy for them to remain straight (their 
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lowest-energy configuration), as the excluded volume 
term would not be a very effective constraint on their 
configuration. In other words, the curvature of any given 
chain for low number densities is expected to be low. As 
a result, it follows that the radius of the sphere on which 
the integration is to be carried out is quite large, com- 
pared to the monomer length. The surface is in this sense 
effectively flat, so that the integration can be carried out 
in a flat two-dimensional manifold. Equation(|^) can now 
be evaluated easily: 
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where b'^ = 3b'^/£p. Now this expression for the propa- 
gator is seen to be analogous to the propagator Gqm for 
the Schroedinger equation in two dimensions, viz.. 
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where r is the time, f is the position vector, and m is the 
particle mass. The mapping between the two models is 
given by TO = -p^, and ir — > —q. 

Since one can use Feynman's formulation to connect 
the quantum-mechanical propagator to a wave function, 
it is natural to ask whether, working in imaginary time, 
one can relate the Green's function in Eqn.(y) to a func- 
tion that plays the role of a probability amplitude. It is 
advantageous to do so, since it permits the formulation of 
a non-trivial model of interacting semi-flexible polymers. 

Our_ formulation uses a standard result from field 
theoryQ to express Eqn.(H) in an alternative form through 
the use of an auxiliary complex field 'ip{t,n), and an as- 
sociated functional Tq: 

goit{N),t{0)) ^ J V^iP rit,N) V'(M) exp(-/3^o) 
(3To - J dnj d'^'Hrii n) (^5„ - ^V^^ ^Pit, n) (5) 

where 9„ = d/dn, 6'^ = 3&^/A, and A — ip/b. 

Even at this non-interacting level the new formulation 
has an advantage over the preceding one in that it is 
a compact way of representing many independent poly- 
mers, since one can easily derive a 2p-point correlation 
functions as a product of p factors of Go by starting with 
the partition function: 



(6) 



Excluded volume effects, representing the simplest 
form of interaction between semi-flexible polymers can 
be expressed by the following extension: 
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where v is the usual excluded volume interaction pa- 
rameter and is a Lagrange multiplier used to con- 
serve the number of monomer segments in the system. 
The model was written down in analogy with the quar- 
tic quantum electrodynamic theory, where an interac- 
tion term local in space and time suffices to describe the 
physics, as long as the coupling constant is chosen cor- 
rectly. It can be thought of as a pseudopotential model 
which is local in the independent variables. 

An Euler-Lagrange equation obtained by extremizing 
(3J- with respect to ip* . An equation of continuity can 
be obtained as in Quantum mechanics from it. \^\'^ can 
now be interpreted as a probability density: 
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While this interpretation of It/i^ as the probability den- 
sity is obtained within a mean field approximation, it 
should not be expected to break down when we consider 
perturbative fluctuations around the mean field. 

The chemical potential can be estimated by consid- 
ering the homogeneous solution of the Euler-Lagrange 
equation, which leads to fi = cqv, where cq is the mean 
number density. 

Now let us re-examine Eqn.(Q). The integrand of /SAT 
is seen to exhibit a double minimum potential as a func- 
tion of ji/'p. As /J, ^ 0, or equivalently cq 0, it can be 
seen that the double minimum vanishes. Thus this 
theory is seen to be analogous to a Ginzburg-Landau the- 
ory used to describe the approach of a thermodynamic 
system to a critical point. Gonsequently, for sufficiently 
small number densities, a system of semi-flexible poly- 
mers approaches a critical point, where fluctuations be- 
come important. We shall account for them by appealing 
to Rcnormalization Group arguments. 

Our interest is,in systems composed of long chains. As 
done previouslyn, we use 9„ —iuj, upon invoking a 
Fourier transform with respect to n. Further, as usual, 
—iuj = \p — q\~^, where \p — q\ is the separation of units 
along the chain. The truncated Hamiltonian that will 
now be considered (for fi —> 0^) is: 
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where N is the chain length in the system. This trun- 
cation allows us to focus on the behavior of semi-flexible 
segments separated by N monomers, e.g., the end-to-end 
distance. Now the entire Wilson formalism for dealing 
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with second order phase transitions can be brought to 
bear on the problem, and it follows that the character- 
istic length scale in the system ^ possesses the following 
scaling form, since d — 1 = 2: 

e - N (10) 

From de Genne's workB and our owni, ^ can be identi- 
fied with the radius of gyration. If the system was com- 
posed of freely extensible chains, then the restriction to 
d~l dimensions would be unnecessary, and we obtain the 
result that Rg ~ ]\[-0-63...^ j|. ^^^^ turn out to be useful 
to study extensions of this theory to non-integer dimen- 
sions, by restriciting the requisite integration to c? — e, 
where 1 > e > indicates varying degrees of extensibil- 
ity of the chains. 

As the number density increases, one can visualize 
these worm-like chains beginning to get entangled, while 
attempting to retain their low-energy linear conforma- 
tion, by weaving themselves into a three-dimensional pile 
of enmeshed semi-rigid needles. To describe this transi- 
tion, it suffices to consider the lowest order vertex cor- 
rection, which arises from the first order bubble diagram. 
It is convenient in this calculation to scale all lengths in 
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the theory by Cq , where cq is the average number con- 
centration of segments in the system: 
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where m — 3/6'^, y = —iio — 2/i, and the 
required frequency integral has been performed by 
the method of residues by closing the contour in 
the upper half-plane, and the denominator, viz., 
{jj + k-^ / 2m + p^ / m + {pk / m) cos 4>) has been exponen- 
tiated using Feynman's trick. The relevant integrals can 
be performed simply, with the result that: 
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Notice that as fc ^ oo, r(fc,a;) — » 0. This is as it 
should be, since F is simply the lowest order density- 
density structure factor. 

In fact, had we done the required frequency in- 
tegral in Eqn.(|ll|) by closing the contour in the 
lower half-plane, we would have obtained, after 
performing dimensional regularization, T{k^Lo) = 
(ll) (V'(l) -ln[m(?/-h3fcV(4TO))]), where ip{l) is a 
digamma function. This form does not vanish for 
large wavenumbers, as a density-density structure factor 
should. 



It follows immediately from Eqn. (|Tl|) and Eqn.([T^) 
that as the number density of segments increases, the 
effective coupling constant decreases. Physically this 
means that the semi-flexible polymers begin to behave 
as if they were increasingly independent. This is a di- 
rect generalization of Flory's theorem for flexible polymer 
melts. 

With Eqn.(p^ in hand, it is possible to define a beta- 
function as is done in Renormalization Group theory, to 
study the scaling properties of the renormalized coupling 
constant ur. 

In the long wave-length limit, the structure factor 
F(fc — 0,cl') diverges logarithmically as y ^ 0+. Note 
also that y = {—iuj — 2/i) = (1/A^ — 2/i) plays the role of 
an effective inverse chain length. We would like to show 
that this is a signal for the onset of entanglement. For 
this purpose we rescale a^, and other relevant quantities, 
for convenience: 
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where y also plays the role of an inverse length scale. It 
has been introduced in order to facilitate a formulation 
of a scaling equation reminiscent of Wilson's Renormal- 
ization Group arguments. A scaling function (3{aR) can 
now be defined: 
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aR = 1 is an infra-red stable fixed point of the transfor- 
mation (in terms of 3^) , and signifies a continuous phase 
transition. This second order phase transition is different 
from the conventional one in that the structure factor di- 
verges merely logarithmically. To connect this transition 
with the onset of entanglement in semi-fiexible polymers, 
we note that the system of semi-flexible polymers can be 
described by the following effective energy functional: 

effective ~ 

d'^-^kduj ^-/ic(fc,w) + ic*(fc,w)r(fc,w)"ic(fc,w)j (15) 

The partition function can be easily computed, 
since only Gaussian path integrals are involved. The 
Helmholtz free energy A follows from this: 



A = -kT\Yi{Z) 
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Using P=- and B ^ (||) , we find: 
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It follows that the bulk modulus exhibits a sharp in- 
crease as 1/N — 2/z ^ O"*", as it must near the onset of 
entanglement. Thus we can identify the phase transition 
located using a Renormalization Group argument with 
the onset of entanglement in a system of semi-flexible 
polymers. The approximations utilized here are insuf- 
ficient to describe the mechanical properties of the sys- 
tem beyond thcr-onset of entanglement. The data pro- 
vided by HinnertI shows evidence of a sharp increase of 
the modulus as entanglement is approached, in a solu- 
tion of actin. There appear fluctuations in the data near 
the transition point. This is qualitatively consistent with 
the (logarithmic) divergence of the structure factor at 
the transition point, which implies large fluctuations in 
the density-density correlation function. The theoreti- 
cal curvature is also consistent with experimental data. 
Note that the number density dependence of the bulk 
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modulus is ~ Cq , similar to the concientration power 
law dependence quoted by Hinner et aI3. It is difficult 
to perform a finer comparison with these data, since the 
experiments were performed for solutions of actin, and 
the theory presented here applies to melts. It should be 
possible to use melts of semi-flexible polymers, where the 
rigidity could stem from steric hindrance, for example, to 
test the predictions of this model. 

Hinner et alEI indicate that in the regime of entan- 
glement, past the critical point discussed above, semi- 
flexible polymers appear to exist in tube-like cages. In 
this regime, away from the critical point just discussed, 
mean-field theory provides an adequate description. To 
investigate the structure of these tubes, let us consider 
the Euler-Lagrange which follows from extremizing J- 
with respect to tp*, written down in dimensionless form, 
in cylindrical geometry: 



ld_ 



dp 



i^s{p? = (18) 



where the radial variable p = t{n) ■ p, with p being the 
unit vector in the radial direction, and s is an integer 
since a functional form exp(is0) was assumed to capture 
the azimuthal dependence of the field. The unit of length 
used for scahng is: 



63/2 



(19) 



For p 0+, ipip) ~ Ap" = A(t{n) ■ p)". Here A is an 
amplitude that can be adjusted in numerical calculations 



to match smoothlyB, over a distance ^ a, the asymptotic 
solution for p ^ oo, Viz., '>p{p) ~ ±\J 1 — / {t{n) ■ pY . 

Physically, what this mean-field solution indicates is 
that due to the inherent stiffness of the polymer, it 
cannot be wound up too tightly. In other words, for 
t{n) ■ p = dnR{n) • p — > O'^, the probability li/jp of finding 
a polymer in that tightly wound position is vanishingly 
small. On the other hand, we do not have free polymers. 
There is an excluded volume interaction term which does 
not allow the semi-flexible polymers to unfurl in a lin- 
ear fashion, but are instead forced to stay coiled, due 
to the presence of other polymers, over a distance ~ a. 
This healing distance a can be regarded as a tube ra- 
dius. Equation (p^ yields an analytic dependence of this 
tube radius on the number concentration of segments, 
persistence length, etc. If £p is of the order of the length 
of polymer chain, then the tube radius is proportional 
to the inverse square-root of the chain length. Since the 
number density of segments cannot exceed 1/v, it follows 
that the smallest value of a attainable is b^'"^/ y^2£p. 

We have used field-theoretic techniques to provide an 
insight into the physics of semi-flexible polymer melts. It 
will be interesting to see how this theory can be extended 
to study polyelectrolytes. 
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